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Abstract. Analyzing (e, e'p) experimental data involves corrections for radiative effects which change the 
interaction kinematics and which have to be carefully considered in order to obtain the desired accuracy. 
Missing momentum and energy due to bremsstrahlung have so far often been incorporated into the simula- 
tions and the experimental analyses using the peaking approximation. It assumes that all bremsstrahlung 
is emitted in the direction of the radiating particle. In this article we introduce a full angular Monte Carlo 
simulation method which overcomes this approximation. As a test, the angular distribution of the brems- 
strahlung photons is reconstructed from H(e, e'p) data. Its width is found to be underestimated by the 
peaking approximation and described much better by the approach developed in this work. The impact 
of the soft-photon approximation on the photon angular distribution is found to be minor as compared to 
the impact of the peaking approximation. 
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1 Introduction 

Much of our knowledge about nuclear structure, e.g. the 
momentum distribution of nucleons in the nucleus, is based 
on (e, e'p) experiments. Currently several such experiments 
are carried out at the Thomas Jefferson National Acceler- 
ator Facility (TJNAF) in Newport News and at the Mainz 
Microtron (MAMI) in Mainz, taking data at high initial 
momenta and removal energies in particular. These exper- 
iments are aiming at a deeper understanding e.g. of short- 
range correlations in nuclei and their results are used to 
check important ingredients of modern many-body theo- 
ries. 



ably; thereby also the momentum transfer between the two 
particles is changed. This phenomenon has been studied in 
a number of papers [^HTOI?^^ 

H7lfT8lfT9] both for inclusive and exclusive electron scat- 
tering experiments. 

Radiative corrections to (e, e'p) scattering, including 
bremsstrahlung, vertex corrections, and vacuum polariza- 
tion (see Feynman diagrams in fig. [1]) can in principle be 
calculated exactly in (pure) QED and to a good accuracy 
also including hadronic loops. But for practical purposes 
several approximations are usually employed when cor- 
recting experimental data for radiative effects [20] . 



In (e, e'p) experiments all particles involved are subject 
to the emission of bremsstrahlung. On the one hand, con- 
sideration of bremsstrahlung contributions is necessary to 
renormalisc the higher-order QED amplitudes. In the sec- 
ond order, divergences from the bremsstrahlung diagrams 
cancel with those resulting from vertex corrections as has 
already been shown by Schwinger in 1949 [T] for electrons 
scattering off an external potential and for electron-proton 
scattering by Tsai in 1961 [2 , including divergent contri- 
butions from the two-photon exchange (TPE) diagrams. 
On the other hand bremsstrahlung modifies the cross sec- 
tion integrated over finite intervals of energy loss. Brems- 
strahlung photons can be so energetic that they influence 
the electron's and the proton's three-momenta consider- 



One of them is the soft-photon approximation (SPA). 
It makes use of the fact that in the limit where ui° — > a 
bremsstrahlung photon with energy uj° has neither a kine- 
matic effect on the scattering process nor an effect on the 
QED propagators and amplitudes. Then the SPA cross 
section factorizes into the elastic first-order Born cross 
section times the probability for emitting a bremsstrah- 
lung photon with vanishing energy. Analysis procedures 
for (e, e'p) experiments make use of the SPA f2U] because 
it simplifies the calculation of multi-photon bremsstrah- 
lung considerably [TS] , 

Multi-photon bremsstrahlung has to be included into 
electron scattering data analysis [TSll2Tj in order to both 
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impose the physical asymptotic behaviour on the cross 
section and to achieve percent level accuracy. While these 
higher-order bremsstrahlung contributions are, in princi- 
ple, also computable exactly in QED, their evaluation has 
to be truncated for practical purposes. The SPA is con- 
venient, because it allows for straight-forward inclusion of 
multi-photon bremsstrahlung into data analysis to all or- 
ders [IS], as we will see in sect. [3] 

The SPA is only valid in the limit of vanishing brems- 
strahlung photon energy; however in radiative correction 
procedures the SPA is applied to photons with energies 
of up to several hundred MeV. The question arises up to 
which photon energies the SPA can be considered a good 
approximation. This paper will not answer that question 
albeit we will present indications that one of the physical 
observables (the missing energy) does exhibit sensitivity 
to shortcomings of the SPA. 

(e, e'p) data analyses do not employ the 'pure' SPA 
(called pSPA in the remainder of this paper) described 
above. Usually the SPA is modified such that it (i) takes 
into account kinematic effects due to emission of finite 
energy bremsstrahlung photons; and it (ii) evaluates the 
form factors at a modified value of the momentum trans- 
fer of the virtual exchanged photon, q. We will refer to 
this 'modified' SPA as the mSPA. In mSPA each particle 
emitting bremsstrahlung is put onto the mass shell. 

The SPA neglects the proton structure at the brems- 
strahlung vertex. But, as has been shown in ref. [15] . for 
the kinematic settings considered in the present manu- 
script, the influence of the proton structure at the brems- 
strahlung vertex is not important. 

The other approximation used in radiative correction 
procedures is the peaking approximation (PA). Most of 
the bremsstrahlung photons from the electron are emit- 
ted either in the direction of the incoming (e) or outgoing 
electron (e') and one can observe two radiation peaks at 
the respective angles. The proton (j>') bremsstrahlung is 
much less peaked. At very high momentum transfers one 
can see a bump (rather than a peak) in its direction, too 
(see fig. [2]). The PA, first proposed for (e, e') experiments 
by L. I. Schiff [3J in 1952, makes use of this observation by 
assuming that all radiation goes either in the direction of 
the incoming electron, or the scattered electron. With the 
advent of coincidence experiments the PA was extended 
to (e, e'p) data [18], assuming that the proton bremsstrah- 
lung was peaked, too. The PA projects the non-peaked 
contributions to the bremsstrahlung photon angular dis- 
tribution onto the three peaks. Especially between the two 
radiation peaks due to electron bremsstrahlung the dis- 
crepancy with data becomes large (see fig. [3]), limiting the 
accuracy of (e, e'p) data analyses [SJE]. 

The purpose of this paper is to remove the PA from 
(e, e'p) data analyses. The need for the removal of the 
PA became evident when looking at the bremsstrahlung 




Fig. 1. The Feynman diagrams contributing to the internal 
radiative corrections together with the Born amplitude (top) 
and the diagrams contributing to the bremsstrahlung ampli- 
tude (below). The label 'iradc' is short hand for internal ra- 
diative corrections and 'brems' stands for bremsstrahlung. 



photon angular distribution in H(e, e'p) experiments (see 
fig. [3]). In this paper, we introduce a full angular Monte 
Carlo (FAMC) method which generates multi-photon brems- 
strahlung events according to the mSPA photon angular 
distribution. A similar FAMC code for (e, e'p) experiments 
has been described in ref. [19]. But it has not been in- 
serted into any data analysis codes nor does it handle 
multi-photon bremsstrahlung. In connection with virtual 
Compton scattering ref. |22j introduces a numerical calcu- 
lation of radiative corrections beyond the PA, but it con- 
siders single-photon emission only, whereas multi-photon 
contributions are large. To check our results against ex- 
perimental data we use the SIMC analysis code [20] for Hall 
C at TJNAF and E97-006 experimental data [23] . 

While we do not want to anticipate the results from 
sec. [5] at this stage we do state here on a preliminary basis 
that removing the PA can only be a first step on the way 
to an improved calculation not relying on the SPA. For 
beam energies envisaged for the TJNAF upgrade, a cal- 
culation going beyond the SPA might become necessary, 
albeit an exact multi-photon bremsstrahlung calculation 
is impracticable. 

This paper is organized as follows: In sect. [2] we in- 
troduce the bremsstrahlung cross section including multi- 
photon bremsstrahlung, discussing the QED divergences. 
Our calculation partially follows ref. [IS], as the result- 
ing equations form the basis for our FAMC calculation. 
In sect. [3] we extend this approach to a FAMC simula- 
tion allowing for any number of bremsstrahlung photons 
emitted into the full solid angle according to the full an- 
gular distribution. In sect. [4] we compare the results of the 
FAMC simulation to the PA using the SIMC code, and 
in sect. Owe discuss scope and validity of the SPA, com- 
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Fig. 2. Single- photon angular distribution of bremsstrahlung 
using the SPA. The full line (1) (black) shows the full angu- 
lar distribution A(O^) from eq. @, the dotted line (2) (red) 
considers pure electron contributions (ee) only, the dashed line 
(3) (green) shows the pure proton contribution (pp). The in- 
terference term is partially negative and is shown in the inset 
graph. It is small compared to the other contributions. The 
momentum transfer is Q 2 = 15 GeV 2 , the entire kinematic set- 
ting can be found in table [1] The dip for the proton is due 
to the fact that a massive particle cannot radiate a photon in 
forward direction. The same is true for electrons but the width 
of the respective dip is extremely narrow 18,22 . 
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15 GeV 2 
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21.00 GeV 
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Table 1. The kinematic setting used in fig. [2] 
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Fig. 3. Bremsstrahlung angular distribution for the H(e, e'p) 
reaction at the kinematics given in table [2] The solid curve (3) 
(black) shows the measured experimental angular distribution 
of bremsstrahlung. The experimental photon angle # 7 is recon- 



structed from the missing momentum according to eq. (|22l) . 
The dotted line (1) (green) and the dashed line (2) (red) cor- 
respond to a Monte Carlo simulation based on the PA and 
take into account detector resolution. The red line takes into 
account the emission of bremsstrahlung from e, e', and p' , 
whereas the green line only allows for bremsstrahlung emit- 
ted solely from either e, e', or p' . The solid line (4) (blue) in 
the vicinity of the proton direction simulates a deficiency of 
the apparatus (proton punch-through effects). The PA over- 
estimates the two photon peaks in the two electron directions 
(e and e') and it underestimates the bremsstrahlung between 
these two peaks as well as between the peaks in incident elec- 
tron direction e and in proton direction p' and the radiation 
on the large-angle tail of the e'-peak. Both data and simula- 
tions account for luminosities and detector efficiencies, and no 
arbitrary normalization factor is needed. 



Q 2 2 GeV 2 





3.120 GeV 


|k'| 


2.050 GeV 


Ip'l 


1.700 GeV 


6>e 


32.5° 


0p 


-38.5° 



Table 2. The kinematic setting at which the standard ra- 
diative corrections and the FAMC simulation are compared to 
data. 



paring it to the exact QED calculation for single-photon 
bremsstrahlung from the electron (which will be called 'I7 
calculation'). 

The speed of light has been set to c = 1 throughout 
the paper. 



2 Bremsstrahlung cross section 

In order to obtain the electron-proton cross section to or- 
der a 2 including bremsstrahlung with energy less than 

^max ; 

< W max ) , (1) 

where oj° is the bremsstrahlung photon energy, the ampli- 
tudes depicted in fig. Q] are considered. The four brems- 
strahlung diagrams contributing to A^broms are divergent 
in the limit of vanishing bremsstrahlung photon energy u>° . 
These divergences cancel the ones both from the TPE dia- 
gramtQ and the vertex corrections [T] . The TPE diagrams 

The divergences from the TPE diagrams cancel with the 
one from the electron-proton bremsstrahlung interference term 
which appears after squaring the full scattering amplitude. 
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are special cases. While consideration of their divergent 
pieces is necessary in order to remove all divergences from 
the scattering amplitudes, their finite contributions are 
known to be negligible in electron scattering experiments 
[251I25U25U57] unless a very small L(T)-contribution is de- 
termined via an LT-separation. Mo and Tsai calculated 
the TPE diagrams approximately using only the nucleon 
intermediate state in the limit where one of the two ex- 
changed photons has zero momentum. They applied this 
approximation both in the numerator and in the denom- 
inator of the fermion propagator [5]. Maximon and Tjon 
improved this calculation by removing this approximation 
from the denominator of the fermion propagator [12pi5j. 
And Blunden, Melnitchouk, and Tjon did the calculation 
using the full propagator [28jl29j ■ According to ref. [28] a 
model-dependent calculation of the influence of the TPE 
yields effects of the order of 1-2% for the kinematic set- 
tings considered in the present paper. Most (e, e'p) anal- 
ysis codes follow the calculation by Mo and Tsai [T5II2T)] . 

The SPA allows us to approximate the four brems- 
strahlung diagrams by a product of the Born amplitude 
times a correction factor. In SPA, e.g. the amplitude for 
incident electron bremsstrahlung can be approximated as 



M c 



eM {1) 

w [ cp 



e ■ k 
to ■ k 



0) 



This amplitude corresponds to the second Feynman dia- 
gram of A^broms in fig. [2 Mep is the first-order Born 
amplitude, u> = u) (1, 1, i7 7 ) the photon four- momentum, 
e is the bremsstrahlung photon helicity vector, and k = 
(fc°,k) the incident electron's four-momentum. The four- 
momentum of the scattered electron will be denoted as 
k' = (fc'°,k'), and for the proton we will use p = (p°,p) 
(incoming) and p' = (p'°,p') (outgoing). 

Evaluating the Feynman diagrams in the SPA one can 
show |18) that the cross section for single-photon brems- 
strahlung is 



da 



dn c df2rydui° dQ c ui° 



(3) 



where fl 1 denotes the bremsstrahlung photon angles and 
da^ /dH e is the Born cross section. In the cross section 
(|3|) the dependences on photon energy and photon angle 
factorize and 



A(J? 7 ) 



P 



P 



4-7T 2 \_ LO ■ k' ' LU ■ p' ' LU ■ k LU ■ p 

does not depend on the photon energy. 



(4) 



Integrating over photon angles and energies the total 
cross section for emitting a photon with energy smaller 
than w max can be written as 



^( W ° 







da^ 

dfir 



"[1 ^brems(^max) ^iradc] • 



(5) 



The necessary integration techniques can be found in [30J , 
the remaining calculations are explicitly carried out in 
ref. |18j . The contributions from vertex correction and vac- 
uum polarization (the internal radiative corrections) are 
included in 



radc 



2a 



where the vacuum polarization contribution is 



i 

3^ 



■ log — 

my 



, (6) 



(7) 



in the ultra- relativistic (UR) limit. This expression does 
not only contain electron-positron loops but also heavier 
lepton and light quark-anti-quark loops, mi denoting their 
respective masses. The bremsstrahlung is contained in 



^brcms (^n 



log 
log 

1 

2 



lie! Ik' 



P°p'° 

max 



log — 
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1 



x log 

1 



Ik' I 



log 



Ikllk'l 



log 



|k||k'| 

M 2 



log 



max 

M 
Ik'i 



log 



Ik'l 



(8) 



also given in the UR limit. The single-photon cross section 
([5| is still divergent in the limit of vanishing a; max . By 
taking into account higher-order bremsstrah-lung (multi- 
photon bremsstrahlung) this divergency is rendered finite 
4.21 and, at the same time, experimental accuracy is 
enhanced. It was first shown in ref. [4] that in fact all 
orders of bremsstrahlung contributions can be considered 
by just exponentiating the bremsstrahlung term in the 
cross section ([5]), yielding 

da da^ 

J^i^i < W max ) = -^-CXp[-<5brcms(w m ax)][l - tirade] • 

(9) 

The index i indicates that an infinite number of photons, 
each with an energy less than w max , is emitted. Exponen- 
tiating <5brcms leads to the correct asymptotic behaviour of 
the cross section (J5j> as w max — > 0. 

We now consider the cross section for emitting n pho- 
tons with an energy larger than an artifically introduced 
energy cut-off parametet E m - la together with multi-photon 
emission of photons with individual energies less than the 
energy cut-off E min [IS], 



da(n, E min ) 



df2 c diOidf2i . . .duj^dfin 



df2 Q 



exp[-(5 br 

^iradc) 



s(-^min)] 
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1 A{S2 1 ) A{Q n ) 



x9(u J ° 1 -E min )...6(u J °-E min ). 

(10) 

In order to integrate this cross section, let us introduce where 
an "acceptance function" x^(!? e ,W] ) , (7\, . . . , u;°, f2 n ), de- 
pending on the kinematical variables of the scattered elec- 
tron and the n photons with energies larger than -E m m- 
Multiplying eq. (fTOf with Xai integrating over all photon 
energies up to an upper boundary E max > E m - m , chosen 
large enough to include all photons, where \\ i s non-zero, 
we obtain the cross section 

1JT ^ = ^ _eX P[ _<5 brem 8 (-Smin)](l ~ tirade) 



oo 



n r E : 



n=0 i=l 



13 



(11) 



This way, cross sections with very general restrictions on, 
e.g., the energy loss by the photons can be calculated. 
The acceptance function xa < 1 is the probability for the 



event to be counted in 



[xa]- 



Proceeding towards a FAMC calculation we evaluate 
the cross section (fTTj) using Monte Carlo integration tech- 
niques. In order to rewrite the integral into the standard 
Monte Carlo sum over randomly selected values, we need 
to rewrite the cross section in terms of probability density 
functions (PDFs) for the photon variables. The PDFs are 
then used to generate the random values. 

As we want to keep the shape of Xa general, it is kept 
in the expression for the cross section and is implemented 
in the Monte Carlo generator by a standard rejection al- 
gorithm. The integral is easily converted into the shape 
demanded by the Monte Carlo generation by multiplying 
and dividing it with its value for Xa = 1- 

The value of the angular integral (without any xa) 



A = J dQyAify) 



(12) 



is independent of the photon energy and A{fi 1 ) is the an- 
gular distribution from eq. @ which is plotted in fig. [2] 
(for = 0) for a sample kinematic configuration (see table 
[1]). The integrals over the lu® can be trivially solved and 
the product over i just yields a power of n. 

Re-writing the cross section in terms of the PDFs for 
the bremsstrahlung photon energies li®, the angular dis- 
tribution and their multiplicity n, leads to 



da 
d£} c 



dn c 

A log 



exp[— £brems(-Eroin)](l - ^iradc) 



xe 



n=0 



PDF(n) 



f[ J du,°PDF(a;°) 



c^PDF^Xa , 



PDF( W °) 



is the PDF for the photon energies wf, 

PDF(/?0 = ^ 
is the one for the angular distribution; and 



(13) 



(14) 



PDF(rc) = — 



A log 



EL 



-A log 



(15) 



(16) 



is the PDF for the photon multiplicity, which is just a 
Poisson distribution. The total cross section (fT5]) does 
not depend on -E m i n because ^brems(-E'min) and A are such 
that E lnin cancels, as has been shown by R. Ent et al. in 
ref. [IS]. 

In the same reference [TB] it is shown that the cross 
section for emitting several photons each with energy less 
than a cut-off E m - m , 



da 
dQ e 



(cj° < -E m ; n ) ■ 



(17) 



is, within a correction of order a 2 , the same for the case 
where instead of the individual photon energies the sum 
of the energies of all bremsstrahlung photons is smaller 
than the cut-off, 



< E„ 



d<T (, ,0 ^ T? 



>)[i + e>(« 2 )] 



(18) 



Therefore the cross section (fl~3)) . within an order a 2 correc- 
tion, can be regarded as the cross section for multi-photon 
emission below a small cut-off i? mm , the sum of these soft 
photons being E m in, along with the emission of n hard 
photons with energies above E m - ln and below £ max . The 
dependence on -E m i n cancels [IS]. For practical purposes, 
in (e, e'p) data analyses, -Emin is often set to a value below 
the detector resolution, and the bremsstrahlung photons 
below that small cut-off are not considered in the anal- 
yses, as even their sum will not affect the measured re- 
sult: Eq. ([15]) ensures that neglect of these photons only 
amounts to missing energies below the detector resolu- 
tion within an order a 2 correction. For the value of E max 
we can always use the total energy of the incoming elec- 
tron, as a reasonable acceptance function disallows events, 
where the photons have together more energy than the to- 
tal energy available. 
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Starting with the cross section in (fl~0|) we can also ob- 
tain differential cross sections. For example, the cross sec- 
tion differential in the total energy of all emitted brems- 
strahlung photons E to t, dn d ^ Et ; , is calculated by choos- 
ing xa = 5 (Er=i <■ 



1500 



^tot 

^tot). yielding 



da 



dQ r dEt. 



dflr. 



CXp[-5 b rems(£ ; min)](l ~ ^iradc) 



xe 



Mfe) ^PDF(n) 



71=0 



(19) 



where we have made use of the fact that in this case the 
integration over the angular variables can be done triv- 
ially. 

In the Monte Carlo simulation events are generated 
according to the PDFs (QUI, (TH]), and |(TBJ) and the results 
are binned in the vicinity of 



(20) 



As discussed above, n, Qi, and uif have to be generated 
according to the PDFs (fT4")) to ([To]). While this approach 
has been used in (e, e'p) data analysis codes (in connection 
with the PA) and is well-established [20], an (e, e'p) data 
analysis code using a multi-photon FAMC simulation is 
novel. In the next section we will describe how the photon 
angular distribution (HJ) is generated. 



3 Full angular Monte Carlo simulation 

For H(e, e'p) data the PA exhibits its limitations espe- 
cially in the middle between the two radiation peaks in 
e and e' direction^ where it underestimates the strength 
of the bremsstrahlung. The same is true for the region 
between the peaks in e and p' directions. In this section 
we introduce a FAMC simulation for multi-photon brems- 
strahlung following the SPA distribution in eq. ([4]) in order 
to see whether this cures the problem. In order to generate 
Monte Carlo events according to the angular distribution 
(01 we need a set of invertible envelope curves Ai(f}~) 
which limit A(f2-y) from above, 



(21) 



for all photon angles Q 1 . In order to obtain the exact dis- 
tribution ^ from the envelope curves we then employ a 



2 Only for H(e, e'p) the missing-momentum vector and thus 
its direction is solely due to emission of bremsstrahlung photons 
below pion threshold. 
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Fig. 4. Angular distribution of bremsstrahlung photons. The 
dotted line (1) (green) and the dashed line (2) (red) repre- 
sent the simulations using the PA (the difference between the 
two curves is explained in the caption to fig. [3} , and the solid 
line (3) (black) shows the measured (reconstructed) experimen- 
tal angular distribution, as in fig. [3] The dot-dashed curve (4) 
(blue) represents our FAMC simulation (not using the PA; em- 
ploying method 3 described in the text; see fig. [BJ. The peaks 
in e and in e' direction generated by the FAMC simulation (4) 
are broader than the ones from the PA, (1) and (2), and repro- 
duce the data more accurately, especially between the peaks 
in e and e' direction. The height of the e' peak is slightly 
underestimated by the FAMC simulation (4), as well as the 
large-angle tail beyond the e'-peak. For details concerning the 
proton direction see fig. O As in fig.|3]both data and simlations 
are corrected for luminosities and detector efficiencies, and no 
arbitrary normalization was introduced. 



standard rejection algorithm. 

The envelope curve chosen by us consists of four terms. 
In order to be able to apply the mSPA we need to assign 
each bremsstrahlung photon to one of the particles. Three 
of the four contributions to the envelope curve can unam- 
biguously be assigned to bremsstrahlung from the incom- 
ing and outgoing electron and the outgoing proton. The 
fourth envelope term takes up the remaining part. It is 
an angle-independent distribution at first but shaped by 
the rejection algorithm into a contribution which is given 
by ee interference. There are several 'coin toss' methods 
to choose whether an event created from the interference 
term is assigned to the incoming or the outgoing electron 
or to both. We employed three different ways of dealing 
with the interference term, leading to slightly different re- 
sults. Together with the Monte Carlo photon energy gen- 
eration (|14p and with the photon multiplicity generation 
(HHJ) each of these three ways of dealing with the inter- 
ference term constitutes a Monte Carlo event generation 
method for the interference term. 
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Fig. 5. Angular distributions of bremsstrahlung photons. This 
figure is a cutout of fig. |4] focussing on the region around the 
proton direction. On the small-angle tail of the proton peak 
the FAMC simulation (4) does not entirely overcome the gap 
between simulation and data. However, the region around the 
proton direction is obscured by punch-throughs. 

1. The interference term (being essentially a function of 
the photon angle 6* 7 ) is split into two parts, the 'left 
part' consisting of events with angles closer to 9 and 
the 'right' part with angles closer to zero. Events closer 
to the incoming electron direction ('right') were counted 
for the incoming electron whereas events closer to the 
outgoing electron ('left') direction were counted for the 
latter one. 

2. In addition to method (1) the energy loss generated 
using (|14|) is randomly split between incident and scat- 
tered electron. 

3. The emitted photon is randomly assigned to either the 
incoming or the outgoing electron. 

For the final comparison between the standard brems- 
strahlung treatment (using PA) and our FAMC simulation 
we used the third method as it fitted the reconstructed 
photon distribution most accurately, as we will see in the 
next section. 

Once a bremsstrahlung event has been assigned either 
to the incident electron, the scattered electron or to the 
struck proton, in mSPA the four- momenta, k, k', and p' 
are replaced by k — > k — u>, k' — » k! — ui, or by p' — > p' — u>, 
respectively, and the momentum transfer q 2 is adjusted 
and inserted into the form factors. u> is the four- vector of 
the bremsstrahlung photon. 

To check the results produced with our Monte Carlo 
routine against experimental data we implemented it into 
the simc code [20J developed for Hall C at TJNAF. We 
used a modified version which was used for the E97-006 



9 y (degrees) 

Fig. 6. Photon angular distribution for the three different 
treatments of the interference term (see sec. [3]). Method (1) 
is represented by the dashed line line (red) and it coincides 
with method (2) (not shown). The dotted line (blue) represents 
method (3), the solid line (black, marked with the black square) 
shows the data. Method (3) is found to reproduce the data most 
accurately and is hence used for the FAMC simulation (4) in 

fig.H 

(e, e'p) experiment 23J. Computation times with and with- 
out the new FAMC simulation were similar. 

4 Results 

To test our approach we chose H(e, e'p) kinematics with a 
beam energy of 3.12 GeV. We usually generated 600,000 
successful Monte Carlo events per run to compare PA and 
FAMC simulation. Figs. [4l - [TOl show results for the kine- 
matic setting given in table [2] 

The photon angles shown in figs. [Hand [4] are obtained 
according to the prescription 

6> 7 = arctan ( — J , (22) 

where p my and p mz are the missing momenta in y and in z 
direction, respectively. Our co-ordinate system is the one 
used by SIMC, described in [3"T] . 

As pointed out in the introduction, the PA underesti- 
mates non-peaked radiation especially between the radia- 
tion peaks in the directions of the incident and the scat- 
tered electron J5[[S]. One can see in fig. 2] that the photon 
angular distribution broadens when employing the FAMC 
simulation. The gap between the experimentally deter- 
mined bremsstrahlung distribution and PA (see fig. [3]) be- 
tween the two radiation peaks in e and e' direction is filled. 
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missing energy [GeV] 

Fig. 7. Logarithmic plot of missing-energy distribution. The 
dotted line (black) was obtained using the PA, the solid line 
(red) shows the results obtained with the FAMC simulation. 
The total radiated energy simulated for the latter case is only 
about 0.3% smaller than the one from the PA. Reasons for this 
deviation are given at the end of sec. [4] 



When calculated with our FAMC method, also the 
peak in the proton direction fits the reconstructed brems- 
strahlung data (see fig. [5]). However, this has to be put 
into perspective as the proton bremsstrahlung is obscured 
by a detector related artefact (punch-through effects) such 
that one cannot make a clear statement on the accuracy 
here. 

For the kinematic setting shown in table [2] the ee inter- 
ference term discussed in the previous section was treated 
with method (3). This led to the best agreement with 
data as can be seen in fig. [6l The other two methods also 
improved the angular distribution of the bremsstrahlung 
but exhibited a slightly larger deviation from the data 
concerning the amplitudes of the e and e' peaks. At the 
kinematic setting shown in table[2]the ee interference term 
accounted for roughly 20% of all events. 

Looking at the missing-energy distribution in fig. [7] 
which includes detector resolution and acceptances, we see 
that the total FAMC yield is 0.3% smaller than predicted 
by the PA, while a calculation not taking into account 
detector resolution and acceptances would yield identical 
results for FAMC and PA. The 0.3% difference is well 
within the systematic uncertainty usually attributed to 
the radiation correction. The missing momenta (see fig.© 
generated by the PA and by our FAMC code do hardly 
differ either, as the missing energy distribution. The mo- 
mentum distributions of electron and proton for the kine- 
matics shown in table [2] are also not changed significantly 




i I i i i I i i i I i i -I r 

0.2 0.4 0.6 

missing momentum [GeV] 

Fig. 8. Logarithmic plot of the reconstructed missing- 
momentum distribution. The missing momentum is almost un- 
altered. The FAMC simulation is represented by the solid line 
(red), the PA by the dotted line (black). 




1.8 1.9 2 2.1 

electron momentum [GeV] 

Fig. 9. Electron momentum distribution. The dotted line 
(black) was obtained using the PA, the solid line (red) shows 
the results using the FAMC simulation. 



by the FAMC calculations, as can be seen in figs. [9l and [POl 

As a further check we also looked at kinematic settings 
with both larger and smaller values of Q 2 , while we let the 
beam energy unaltered. We compared again the FAMC 
simulation with the standard radiation code. Looking at 
the total yield in the acceptance we found differences of 
up to 3.0%, the yield of the FAMC simulation usually be- 
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Fig. 10. Proton momentum distribution. FAMC simulation 
(solid line) (red) and peaking calculation (dotted line) (black) 
almost coincide. 



Q7GeV 


0.61* 


1.00* 


2.00 


3.00* 


4.00* 


Ip'j/GeV 
k' /GeV 
A yield 


0.852 

2.74 

+2.5% 


1.13 
2.59 
+0.4% 


1.70 
2.05 
-0.3% 


2.36 
1.52 

-1.5% 


2.92 
0.99 
-3.0% 



Table 3. The differences in the yield between the stan- 
dard analysis code and the FAMC simulation, integrated up 
to 0.7 GeV. The incident electron's energy is k° = 3.120 GeV 
for all kinematic settings. Reasons for the differences are given 
at the end of sec. [4] Missing energy plots for the kinematics 
settings marked with a star are shown in fig. 1111 



ing smaller than the standard analysis yield when going 
to higher momentum transfers and larger for small values 
of Q 2 , as can be seen in table [3] and in fig. [TTJ 

The differences in the total yield shown in table [3j in 
figs. 171 and II II are related to the inappropriate application 
of the SPA. It only shows up when including detector sim- 
ulations into the data analysis. Our FAMC approach is 
more sensitive to problems caused by the SPA than the 
PA at certain kinematic settings. It uncovers a problem 
of the SPA which is suppressed by the PA. Including the 
full angular dependence of bremsstrahlung photons (other 
than the trivial angular dependence of the PA) can some- 
times lead to energy gains for both electron and proton as 
the particles are assumed to be on-shell. Such un-physical 
events are rejected by our code, because they are artefacts 
of the mSPA which corrects for the energy losses due to 
photon emission and assumes on-shell vertices. At some 
kinematic settings the un-physical events described above 
account for a significant fraction of all events, changing 
the total yield. 




-0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



missing energy [GeV] 

Fig. 11. Logarithmic plots of missing-energy distributions 
for the remaining kinematic settings shown in table [3] The 
distribution for Q 2 — 2.00 GeV 2 has already been shown in 
fig. [7] The solid curves (red) represent the FAMC calculation, 
the dotted curves (black) show the PA calculation. 

The PA does not have this particular problem since it 
can fulfil both energy and momentum conservation at the 
same time when assuming massless on-shell electrons. En- 
ergy gains through emission of radiation are not possible. 
The recoiling protons cannot be assumed to be massless, 
of course. But as they only account for a small fraction 
of high energy bremsstrahlung events, they do not change 
the total yield much, neither in the case of the PA nor for 
the FAMC simulation. 



5 The applicability of the SPA 

As we have shown in sect. [2] the SPA simplifies the multi- 
photon bremsstrahlung treatment considerably. In order 
to evaluate its applicability for the kinematic settings con- 
sidered in this paper we now test the SPA by comparing it 
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Fig. 12. Photon energy distribution for the kinematic setting 
shown in table[2] The dotted curve (1) (black) represents pSPA 
where the kinematics used are purely elastic. The dashed curve 
(2) (green) shows the mSPA calculation which has been cor- 
rected for the kinematic changes due to bremsstrahlung. This 
mSPA has been used in the present paper and also in refs. [181 
120] . The solid curve (3) (red) depicts the exact I7 calculation 
for the photon energy. The inset graph shows the difference 
between the calculations in percent. 



0.06 



0.04 S 



T3 




0.02 



0.00 

10 

e y [degrees] 

Fig. 13. Bremsstrahlung angular distribution for single- 
photon emission at the kinematic setting shown in tabled The 
dotted curve (1) (black) represents pSPA, the dashed curve (2) 
(green) shows the mSPA calculation, the solid curve (3) (red) 
depicts the exact I7 calculation for the photon angles. The 
inset graph shows the deviation between the different calcula- 
tions in percent. 



to the exact I7 calculation (omitting proton bremsstrah- 
lung). The integration over the bremsstrahlung photons 
is carried out with our FAMC generator. This renders an 
analytic evaluation of phase space integrals unnecessary. 

Let us first describe how we construct a Monte Carlo 
generator for the exact I7 bremsstrahlung calculation from 
the mSPA Monte Carlo generator described before. Con- 



sider the SPA bremsstrahlung cross section 



OSPA 



\Mg>\*A(w) 

\MV>\*A(u,), 



2w '"' cp 
u doj°da 



(23) 



where we have absorbed the photon energy u>° into the 
SPA angular distribution, 



A(co) 



(24) 



Our FAMC code generates bremsstrahlung events accord- 
ing to the distribution 



-AM. 



Evaluating the phase space integral in Eq. 
FAMC generator we obtain 



(25) 
using our 



CSPA 



events 



U! 



2 A(u)< 



where N is the number of events, and Aiip is the elastic 
first-order Born matrix element. The exact I7 calculation 
(not using the SPA) yields 



(26) 



u>°dLu a df2~, ... ,9 
o ^M^l 2 . 



where 



Ml-y = M C 



M 



cf 



(27) 



(28) 



is the exact QED single-photon electron bremsstrahlung 
amplitude. The cross section (|2"7|) becomes 



(7 1- 



- V 

N ^ 

events 

N ^ 



\M 



11 



1 



2 



\M 



It 



A(u) 



(29) 



in the Monte Carlo formalism. In order to measure the 
applicability of the SPA we assign a weight w cx to each 
event, defined by the ratio of the squared matrix elements 
from cqs. ([26| and ([29]) . re- weighting our FAMC generator 
(SPA), 



\M 



1"! 



\M^\ 2 A{uj) 



We compare that weight with the mSPA weight 



W mo d 



_ |Mi 1 p ) ' mod | 2 J 4 mod ' 



\M$\*A(u) 



(30) 



(31) 
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where 'mod' indicates that these matrix elements in the 
numerator have been calculated for modified kinematics, 
i.e. in mSPA. Finally we define the trivial weight 



w tr 



\M^\ 2 A{lo) 



(32) 



which represents the pSPA, i.e. the pure SPA calculation 
not taking into account any kinematic changes imposed by 
emission of bremsstrahhmg. This last option is the rough- 
est approximation. A version of the mSPA represented by 
weight ([3"T|) is also used in ref . [TH] , combined with the PA. 

Plotting the photon energy distribution (see fig. [T2"]) us- 
ing the three weights (|30|) , ([3"Tf . and ([32]) for each Monte 
Carlo event, we see that the deviation between the mSPA 
calculation and the I7 calculation for the photon energy 
is 4.1% for uj° — 100 MeV. The deviation becomes much 
larger for higher energies, going up to 90% for photon en- 
ergies of uj° = 1000 MeV, the mSPA calculation [1511213] 
overestimating the radiative tail; and one can find even 
larger deviations for different kinematic settings. How- 
ever, bremsstrahlung events with photon energies of se- 
veral hundred MeV are unimportant for the data analyses 
since the particle detectors do not see them. Their mo- 
mentum acceptances usually are limited to ±10% of the 
central (elastic) momenta, depending on what one is look- 
ing for. 

The dotted curve (1) in fig. [T2l shows the pSPA which 
was obtained using weight (|32[) . Most data analysis codes 
make use of some version of the mSPA. The mSPA is closer 
to the exact calculation than pSPA, so it constitutes an 
improvement over pSPA. This finding is in agreement with 
ref. [IE]. 

For our purposes the influence of the SPA on the pho- 
ton angular distribution is more important than the devi- 
ation in the missing energy calculation. Fig. 1131 shows the 
photon angular distribution. The largest deviations occur 
in the vicinity of the peak due to radiation from the inci- 
dent electron, at small values of the angle 6* 7 . While the 
I7 calculation of the bremsstrahlung cross section is sym- 
metric in e and e', the mSPA data analysis procedures are 
not, resulting in asymmetric deviations from the I7 result. 
This can be understood from the energy loss of the inci- 
dent electron, leading to smaller bremsstrahlung energies 
coming from the scattered electron. From fig. [3] we know 
that one critical domain of large discrepancies between 
data and standard simulations using the PA is the region 
in the middle between the two radiation peaks where the 
PA angular distribution falls below the measured distri- 
bution by a factor 2. Fig. [TBI shows that the FAMC cal- 
culation (using mSPA) overestimates the photon angular 
distribution in this region by 13%. Fig. 2] however sug- 
gests, that the FAMC (using mSPA) reproduces the data 
well, especially in the region in the middle between the two 
radiation peaks. Comparing fig. Q] (which includes inter- 
nal and external bremsstrahlung, multi-photon emission, 



finite detector resolution and acceptances, multiple scat- 
tering and other energy losses) and fig. [13] (internal single- 
photon bremsstrahlung only) in the critical region between 
the two electron radiation peaks, we can conclude that the 
SPA impact on the angular distribution is smeared out by 
other sources of inaccuracies like finite detector resolution, 
finite momentum acceptances, and multiple scattering. 

In fig. 0] we saw that the height of the e' peak is slightly 
underestimated by our FAMC simulation, as well as the 
large-angle tail beyond the e' peak. The results shown in 
fig. [Tj2 indicate small differences between the mSPA and 
the exact I7 calculation in the vicinity of the e' peak only. 
And also the difference in fig. [13] at the large-angle tail 
beyond the e! peak is small. At this stage it is not entirely 
clear whether removal of the SPA would affect the photon 
angular distribution in these regions. 



6 Conclusion and Outlook 



Using a FAMC bremsstrahlung calculation at almost no 
extra computational expense improves the treatment of in- 
ternal bremsstrahlung in (e,e'p) experiments. One short- 
coming of the PA, the underestimation of bremsstrahlung 
between the radiation peaks, is solved by our approach. 
We have also shown how the PA can be removed. 

Figs. [7] and [TTJ and table [3] indicated that the FAMC 
simulation exposes problems due to the SPA which are 
hidden when using the PA. And while for the photon an- 
gular distribution the PA may be the dominant source of 
error, we have shown in the previous section that the SPA 
seems to have a sizeable influence on the missing energy 
distribution. 

These two problems with the SPA indicate that it 
would be desirable to also remove the SPA from data anal- 
ysis codes. But there are several problems which have to be 
tackled in order to achieve such an improved calculation. 
An exact multi-photon calculation would be impracticable 
since one would have to insert the QED cross sections into 
the analysis codes for multi-photon bremsstrahlung up to 
arbitrarily high orders. One could instead try to combine 
exact single-photon bremsstrahlung for the hardest brems- 
strahlung photon with SPA multi-photon emission for the 
softer photons in order to improve the present bremsstrah- 
lung treatment. Yet, inclusion of proton bremsstrahlung, 
as in the present paper, seems not to be feasible for such 
calculations. 
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